liep-th/9904116 



April 1999 



o 



X 



Breakdown of Cluster Decomposition in 
Instanton Calculations of the Gluino Condensate 



Timothy J. Hollowood^'^, Valentin V. Khoze^, Weonjong Lee^ 

AND Michael P. Mattis^ 



CN I ^ Theoretical Division T-8, Los Alamos National Laboratory, 

(N ' Los Alamos, NM 87545, USA 

^ '^Department of Physics, University of Durham, 



Durham, DEI 3LE, UK 

^Department of Physics, University of Wales Swansea, 
Swansea, SA2 8PP, UK 



ABSTRACT 



A longstanding puzzle concerns the calculation of the gluino condensate (||^) = cA^ in 
Af = I supersymmetric SU{N) gauge theory: so-called weak-coupling instanton (WCI) cal- 
culations give c = 1, whereas strong-coupling instanton (SCI) calculations give, instead, 
c = 2((A^ — 1)!(3A^— 1)) . By examining correlators of this condensate in arbitrary multi- 
instanton sectors, we cast serious doubt on the SCI calculation of (||^) by showing that an 
essential step — namely cluster decomposition — is invalid. We also show that the addition of 
a so-called Kovner-Shifman vacuum (in which (y|~^) = 0) cannot straightforwardly resolve 
this mismatch. 



I Introduction 



Almost dating back to the development of QCD itself, supersymmetric versions of QCD 
have been closely studied, as tractable laboratories for extracting exact analytic information 
about both perturbative and non-perturbative phenomena in nonabelian gauge theories. One 
outstanding puzzle, unresolved since the mid-1980's, concerns the calculation of the gluino 
condensate (||^) in these models. This is an interesting quantity, as it is a measure of chiral 
symmetry breakdown. In pure supersymmetric Yang-Mills (SYM) theory, by dimensional 
analysis, one expects 

/trA^ 



\l67r 



cA' , (l.l) 



where A is the dynamical scale in the theory (developed by dimensional transmutation as in 
QCD), while c is a numerical constant. Remarkably, there are two approaches in the literature 
for calculating (y|^), each purporting to be exact (i.e., nonrenormalized), but which differ 
in their predictions of the constant c. This disagreement is especially vexing in light of the 
fact that both involve the use of supersymmetric instantons. The first approach, generally 
known as "strong-coupling instanton" (SCI) calculations, was developed in Refs. while 
the second approach, generally known as "weak-coupling instanton" (WCI) calculations, was 
developed in Refs. for self-containedness, both will be reviewed below. 

In this paper, we re-examine this old controversy, using our recently developed methods 
for studying supersymmetric multi-mstantons ||9|-[12|. In particular, by looking at n-point 



correlators (-^^y^^^ ' ' ' ^^1^^) gluino condensate, we will be able to probe arbitrary 

topological numbers k. In a nutshell, our results cast serious doubt on the validity of the SCI 

calculations of the condensate. Specifically, we will demonstrate that an essential technical 

step in the SCI approach, namely the use of cluster decomposition, is invalid. The important 

implications of this observation are as follows. Since cluster decomposition is an essential 

requirement of quantum field theories (with very mild assumptions that are certainly met 

by SYM theory), the exact quantum correlators must have this property. That cluster is 

violated by the instanton-saturated SCI correlators then means that (contrary to claims in 

the literature) the SCI approximation is only giving part of the full answer. Since the SCI 

correlators obey supersymmetric perturbative nonrenormalization theorems ^ , it necessarily 

follows that additional non-perturbative objects must be contributing to the correlators. A 

fuller discussion of this point is given in Sec. VII below; however, categorizing the nature 

of these non-perturbative configurations is beyond the scope of the present paper.0 We 

^See however Ref. where, in a compactified version of the present theory, the important role played 
by monopoles is emphasized. 
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should add that we beheve that, in contrast, the WCI correlators are consistent with cluster 
decomposition. 



In addition, we will address an ingenious, if controversial,^ hypothesis of Shifman's, in 
which the numerical disagreement between the SCI and WCI results is taken as circumstantial 
evidence for the existence of an extra disconnected vacuum in SYM theory in which chiral 
symmetry is unbroken |16, 1^. While this so-called "Kovner-Shifman (KS) vacuum" can 



indeed potentially resolve the disagreement at the 1-instanton sector {k = 1), we will show 
that it fails to do so for the topological sectors with /c > 1. In other words, positing a KS 
vacuum cannot by itself restore the cluster property to the SCI correlators. This discouraging 
finding might be viewed as removing some of the impetus for positing such a vacuum in the 
first place. 

Finally we will present a novel calculation of {^^^^) which relates the Af = I supersym- 
metric models discussed herein to the exactly soluble Seiberg-Witten models with M = 2 
supersymmetry. This calculation is of potential pedagogical interest because it bypasses the 
explicit use of instantons, and instead relies on functional methods. Not surprisingly, it 
recaptures the WCI answer. 

Let us sketch in broad strokes the main differences between the SCI and WCI calculations 
(a more detailed review will follow). For M = 1 supersymmetric SU{N) gauge theory with 
no matter, the leading coefficient of the /9-function is bo = 3N, so that goes like an 
"jV*^ root" of an instanton: A'^ oc exp(— 87r^/(7^A^). This means that a naive 1-instanton 
calculation of (y|^) — in which A is simply replaced by its "classical value" as an adjoint 
fermion zero mode in the instanton background, and all the instanton collective coordinates, 
both bosonic and fermionic, are integrated over — fails; specifically it gives a zero answer, due 
to unsaturated Grassmann integrations. In order to perform a sensible 1-instanton calculation 
of (H^)) two alternative, and necessarily more elaborate, approaches suggest themselves. In 
the SCI approach, one calculates the A^-point correlator of this condensate, which scales 
like exp(— 87r^/(7^), and is indeed nonzero at the 1-instanton level. Furthermore, by a Ward 
identity reviewed in the Appendix, it is independent of the A^ space-time insertion points Xj. 
After performing the requisite collective coordinate integration, one finds: 

trA^(a:i) trA^\ ^ 2^ 
167r2 167r2 / (A^ - 1)! (SA^ - 1) ' ^ ' ' 

In order to extract {^^) from the correlator ( |1.2| ), one then invokes cluster decomposition: 
taking \xi—Xj \ ^ fi^^ where fi is the mass gap in this theory, and remembering the constancy 



2See Ref. and the rebuttal Ref. |T|. 
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(WCI result) . (1.4) 



of the correlator, one replaces the left-hand side of Eq. (|1.2|) simply by {^^)^ ■ The net result 
thus reads: 

where u = 0, . . . , N — 1 indexes the N vacua \u) of the SU{N) theory, and reflects the 
ambiguity in taking the N^^ root of unity. In retrospect (as argued in Refs. 0,§1), the reason 
why the naive calculation of gives zero is that these vacua are being averaged over 

and the phases cancel. 

In contrast, in the WCI approach, one modifies the pure gauge theory by adding matter 
superfields in such a way that (y|~t) itself (rather than a higher-point function thereof) 
receives a nonzero contribution at the 1-instanton level. Next, one decouples these extraneous 
matter fields by giving them a mass M, and taking the joint limit M —* oo and A ^ in 
the manner dictated by renormalization group (RG) decoupling. Matching onto the effective 
low-energy theory without matter gives: 

/ trA^ \ 
\ 167r2 / 

Note that the RG decoupling procedure forces the low-energy theory into one of the A^ 
degenerate vacua \u), which by convention we take to be the one with real phase. The 
nomenclature "strong coupling" versus "weak coupling" used to designate these differing 
approaches refers to the fact that, in the former, as in QCD, the only scale in the problem 
is the dynamical scale A, whereas in the latter, the existence of VEVs Vj of the matter 
superfields permit a standard semiclassical expansion when the dimensionless ratios A/vj are 
all small. (The holomorphic properties of SYM theory then permit the analytic continuation 
of the answer beyond this regime.) 

As mentioned above, it is possible to reconcile the two calculations (|1.3| ) and ( |1.4| ) by 
positing the existence of an extra vacuum IS") in which the condensate vanishes |16|. Specif- 
ically, if p and 1 — p represent the probability weights in the vacuum sector of the theory for 
the standard vacua {!«)}, and for \S), respectively, and if one takes 

^ ^ {N-1)\{3N-1) ' 

then both the 1-instanton results can be understood. Unfortunately, the multi-instanton cal- 
culations presented below show that the mismatch between the SCI and the WCI calculations 
becomes more severe for higher topological number k, and apparently cannot be reconciled 
in this way for k > 1. 
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This paper is organized as follows. In Sees. II and III, respectively, we review the SCI 
and WCI calculations of {^^) for general gauge group SU{N). Also in Sec. Ill we present 
an alternate, non-instanton-based derivation of this condensate, specific to the gauge group 
SU{2), which starts from the Seiberg-Witten solution of the M = 2 model and flows to 
the N = 1 model, recapturing the WCI result. In Sec. IV we discuss cluster decomposition in 
more depth, and motivate Shifman's proposal for reconciling the SCI and WCI calculations 
by postulating an extra vacuum state. Our principal results are described in Sees. V and 
VI, in which (extending the SCI approach) we calculate higher-point functions of the con- 
densate, in the topological sectors A; > 1. In Sec. V we calculate, analytically, the (A;A^)-point 
functions {^^^^^ ■ ■ ■ *''^i6^2^"* ) SU{N) gauge theory for arbitrary instanton number k, but 
to leading order in while in Sec. VI we calculate, numerically, the 4-point function 

( ^'^le^a^^ ^'^lei^^^ ^'"le!-^^'* ^''lel'^^'* ) 2-instanton level for gauge group SU{2). In either case 

our SCI calculations explicitly contradict the hypothesis of cluster decomposition — both with 
and without an extra KS vacuum.^ Concluding comments are made in Sec. VII. 



II Review of the Strong-Coupling Instanton Calcula- 
tion 

Let us review the SCI result for {|trA^(xi) ■ ■ ■ trA^(xAr)), for pure Af = 1 SU{N) gauge 
theory. The calculation for done originally for the SU (2) theory in and then extended to 
the SU{N) theories in P| (see also the very comprehensive review articles H,|T7|). 

The correlator in question is saturated at the 1-instanton level. The gauge-invariant 
collective coordinate integration measure is a suitable generalization of the Bernard measure 
r3 to an A/" = 1 theory, and readsfl 



' {N - 1)\ {N - 2)\ 



j d^a' dp" (p^)2^-^ d^M' d\ d^-^u d^-^u . (2.1) 



Here a'^ is (minus) the 4-position of the instanton and p is its scale size, the Grassmann spinors 

Ai'^ and (a parametrize the supersymmetric and superconformal modes, respectively, of the 

■^The numerical calculation is based on a Monte Carlo integration, which (with our present statistics) is 
incompatible with the clustering result at the 5^ sigma level, and incompatible with the modified clustering 
result due to the incorporation of a KS vacuum (tuned to reconcile the WCI and SCI 1-instanton results), 

at the 11^ sigma level. 

Our choice of notation is dictated by the fc-instanton generalization of this measure, Eq. (|l8|) below. 
Following Rcf. [ pO[ , we correct a factor of two mistake in the normalization of adjoint fcrmion zero modes that 
pervades much of the literature (e.g., Refs. [^D). Hence our final result for the A^-point function, Eq. (1.2), 
differs by 2^ from these references. 
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gluino, and the Grassmann parameters and Uu', u' = 1, . . . ,N — 2, are the superpartners to 
the iso-orientation modes which sweep the instanton through SU{2) subgroups of the SU (N) 
gauge group (note that each z/^/ and Uu' is a Grassmann number rather than a Grassmann 
spinor). The measure includes the Lambda parameter of the Pauh-Villars (PV) scheme which 
at the two-loop level is ^ 

A = ^^(/i)-2/3e-8-V(3^5(^)');, . (2.2) 

Into this measure one inserts Hill where A"(a;) is the most general classical adjoint 

fermion zero mode in the 1-instanton background. In terms of these bosonic and fermionic 
collective coordinates, one derives (see Eq. ( |5.21| ) below): 



o 4 



2^2 



^p^v Zitx {p +y 



(2.3) 



where 



Now let us carry out the Grassmann integrations in Eq. (2J.). Obviously the C, and 
AA.' Grassmann integrations will be saturated from the condensates inserted at two points 
{xj, Xj} chosen from among the insertions Xi, ... , x^v- For each such pair there are three 
contributions to these integrals: 

^ ^l%^f^ {C at x„ M'' at X,) , (2.5a) 

[p^ + yt)[p +yj) 

16 (p2 + y2)2^p2 + ^2)2 

,,2 ,\lyJ{ .,2^2 X at x„ C X at xj) . (2.5c) 



16 ^ ^(p2+ ^2)2(^2 + ^ 



Adding these three contributions gives the simpler expression 

— 36p* (xj — Xj)^ 



(2.6) 



Now we take advantage of the fact that this A^-point function is independent of the Xj (see 
the Appendix), to choose these insertion points for maximum simplicity of the algebra. The 
simplest conceivable such choice, Xj = for all i, turns out to give an ill-defined answer of the 
form "0 X oo" (the zero coming from the Grassmann integrations as follows from Eq. (|2.6| ), 
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and the infinity from divergences in tlie integration due to coincident poles). In order to 
sidestep tliis ambiguity, one cliooses instead: 



xi = ■ ■ ■ = xn-1 = , xn = X . (2.7) 

Tliis clioice is tlie simplest one wliicli gives a well-defined answer with no "0 x oo" ambiguity. 
More ambitiously, one can still perform the calculation even if all the insertion points are 
taken to be arbitrary 0,0; however, we find it convenient for later to take the minimal 
resolution provided by ([2^^ ). From the (xj — XjY dependence in Eq. ( p.6|) , it follows that 
the pair of insertions {xi, Xj} responsible for the {C, A^'} integrations must include the point 
xm = x\ there are — 1 possible such pairs, giving 

-36(iV- l)p8a;2 



(p2 + (x + a')2) (p2 + a'2)^ 



(2.8) 



for these contributions. The remaining Grassmann integrations over {z/, z/} are saturated at 



Xj = 0, and give 



4n2 X N-2 

(^-2)!(7-^^) • (2.9) 



(p2 + a 



Combining the denominators in Eqs. ( p.8|) -( p^ with a Feynman parameter a, 

1 1 (SAT+l)! , a\l-af^-^ 

— da 



(p2 + a'2)3A^-2 ^p2 + + ^,)2)4 3, (3^ _ 3), ^ ^ ^ ^ _ ^)^2)3A^+2 

(2.10) 

and performing the d'^a' integration then yields: 



2 / „2\2N-A 



2iN+2 ^2N-2j^3N rl 

(trA^(xO---trA2(x^)> = _ _ da dp' (p 



X (36(iV- l)pV)(Ar-2)!(4p2)^^2 



3! (3iV - 3)! 3iv(3iV + l)(p2 + a{l - a)^^)'^ (2.11) 

da a'il — a) 



3(3iV - 2) 2^^-! vr^^A^^ 



{N-2)\ 7o 

257V ^2N^3N 



{N -1)\{3N -1) 
in agreement with Eqs. (p..2[)-(P!75D. 
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Ill Review of the Weak- Coupling Instant on Calcula- 
tion 



Next, let us review the WCI calculation of the gluino condensate. As mentioned above, the 
general WCI strategy is to extend the pure gauge theory to include matter content, in such 
a way that (y|~3") receives a nonzero contribution at the 1-instanton level. Decoupling the 
extraneous matter and matching to the low-energy pure gauge theory is then accomplished 
using standard RG prescriptions. Since the precise nature of this extraneous matter is rather 
arbitrary, the WCI calculation really stands for a family of related calculations sharing this 
basic approach, all of which give the same result ( |1.4| ). Calculations of this type were done 



in l5H8i and reviewed in 17 



We will find it efficient to exploit the functional identity (see for example |2l||): 



;trt^> . -8«(Aw.A = ^(aAvv.,) . (3,i; 



Here Weg is the effective superpotential, 

47rz 9 

is the usual complexified coupling, and 



^ = — + (3-2) 



is the RG-invariant 1-loop dynamical scale of the theory. This result comes from writing the 
microscopic gauge theory as 

C = — Im(r [ ^etrWW^) , (3.4) 
An ^ J 

where W"' is the gauge field-strength chiral superfield, and promoting r to a "spurion super- 
field", 

r ^ T{y,9) = T{y) + V2 9''xl{y) + 0'F^{y) . (3.5) 



From Eqs. ( p.4|) -( p3D it trivially follows that 

Sir 



trX') = ^-J—Z , (3.6) 

/ -7 XPrf^\ T{y,9)=T ^ ^ 



where 

Z = I VWe' f'^'''^ (3.7) 
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is the partition function of the microscopic theory, in the generahzed background field 
In order to derive Eq. (|3.1| ) from Eq. (|3.6|) , one assumes that the functional differentiation 
indicated in Eq. ( |3.6| ) formally commutes with the integrating-out of the microscopic degrees 
of freedom. In other words, Z can be re-expressed in terms of the relevant effective chiral 
superfields $i (whatever these may beQ): 

trA^> = 1^7^ ^cfT , (3.8) 

' '7^KT?T(nn\ T(y,9)=T 



where 



Zeff = j p$^e-i''^'^/'^'^^«=«(*-^) . (3.9) 
Equation ( p.l| ) then follows from the observation that dWes/dF'^ = 9^ dWes/dr. 



We now need an explicit expression for the effective superpotential. Following Affleck, 
Dine and Seiberg (ADS) 0, it is convenient to start from SU{N) gauge theory where the 
number of flavors A'^^' is fixed to Np = N — 1. A 1-instanton calculation of the superpotential 
then gives: 

where the flavor indices f,f' = 1, . . . , Np run over the quark superfields. The coefficient of 
the /3-function is, for general A^ and A^, 

bo = 3N-Nf . (3.11) 

The normalization constant for the specific case Np = N — 1 was fixed by an explicit 1- 
instanton calculation, and is simply p0|,p2| Cads = 1- By decoupling the quark flavors one 
at a time, this 1-instanton expression flows into models with Np < N — 1 for which the 
superpotential is no longer a 1-instanton phenomenon. In this way one generalizes Eq. ( p.lO[ ) 



to (see e.g., Refs. ||21 



where 



^N.,N ^ c^^F,^ I ^t^^ I (AT^ < AT - 1) , (3.12) 

detNp[QfQf') 



C^F^^ = N-Np . (3.13) 



In the example culminating in Eq. ( |3.16 ) below, we will find that there are in fact no residual chiral 



superfields so that simply Woff — yVcs{T), whereas in the Seiberg- Witten example ( ^3.17 ) below, the $i 
are the monopole superfields M, M as well as the dual Higgs Ad. 
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Starting from this more general superpotential, let us decouple the remaining quarks, by 
giving them a common VEV v. Viewing Q as an Np x matrix, one assumes: 



/v 

1° 



■■■ o\ 

V ■■■ Oy 






Vo 



V 



0/ 



(3.14) 



The D-flatness condition together with a global gauge rotation gives v = v. Taking |v| — > oo 
then decouples the quarks as well as a subset of the gauge fields, leaving a pure SU {N') gauge 
theory with N' = N — Np and 6o = 3A^'. The 1-loop RG matching prescription reads [pO| : 



A 



Nf,N 



V 



■iN-Np 



Ao. 



N' 



3N' 



(3.15) 



Inputting Eqs. ([CTp- (|37[5|) into Eq. (^J2|) gives: 



cff 



x3N-Nf \ N-Np 



(3.16) 



The desired result (|1.4| ) then follows from Eq. (|3.1| ). 



Note that the starting-point for this WCI calculation, Eq. ( |3.1CI| ), is a bona fide 1-instanton 
calculation. The remaining steps towards the answer involve well-studied path-integral and 
renormalization group manipulations (principally Eq. (p.l|), and Eqs. ( |3.12| )-( |3rT5| ), respec- 
tively). Alternatively, starting again from the functional identity ( p.l|) , we can rederive the 
WCI result ( [1.4|) without any reference to an instanton calculation. Instead, one starts from 
the Seiberg-Witten solution of the Af = 2 model,Q in the presence of a mass deformation 
which breaks the supersymmetry down to A/" = 1. In the strong-coupling domain, in the 
vicinity of the monopole singularity, the superpotential looks like: 



sw 



V2MAdM + mU{A 



D 



(3.17) 



Here the chiral superfields {M , M} describe the monopole multiplet, Ap, is the dual Higgs, 
U is the quantum modulus of the theory (here, in strong coupling, expressed in terms oi Ap 
rather than A), and m is the mass parameter. The F-flatness condition for the vacuum reads 

9Wsw 5Wsw 5>Vsw 







dM 



dM 



dA 



D 



(3.18) 



^For the remainder of the section, we focus on SU{2) gauge theory, and quote weU-known formulae from 
Seiberg and Witten 
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which is solved by 

aj,^{Ao) = 0, (M) = (M) = ( - -^f/'(0)) . (3.19) 
In the vicinity of this solution, the relationship between and m = (f/) is given by 

an = ^ rdx^^= (3.20) 

from which it follows that 



1 V^^^ 



u = Al^- 2tAs^aD + 0{al) . (3.21) 

Here the Seiberg-Witten dynamical scale Ag^ is related to the conventional PV/DR scale 
Ajv'=2 via 



Asw = V2Aj^=2 . (3.22) 

Note that the series ( p.21| ) is not an instanton expansion (i.e., an expansion in Ag^); instantons 
emerge only in the weak-coupling regime, where u is expanded in terms of a = (A) rather 
than an- 



Applying the identity (|3J]) to Wsw using Eqs. ( |3.21| )- (|3^ ) gives the gluino condensate 
in the vacuum (|3.19|) : 

(trA^) = IGTr^mA^^^ . (3.23) 

Next we decouple the adjoint Higgs superfield, by taking m — > oo. In this way we flow to 
the pure A/" = 1 supersymmetric SU (2) gauge theory. The RG matching condition between 
the scale A of the J\f = 1 theory and the scale Aj^^a of the mass-deformed M = 2 theory 
reads |2§: 

m'K^^, = A' . (3.24) 
Substituting Eq. (|3.24| ) into Eq. (|3.23|) once again gives the WCI answer ( |1.4|) . 



IV Comments on Cluster Decomposition 

In this section, we examine the issue of cluster decomposition in the context of the gluino 
condensate. This issue of cluster decomposition is fundamental to a quantum field theory. 
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The clustering property requires that for sufficiently large separations | I, compared 

with the inverse mass gap£| 

{(pi{Xi) ■ ■ ■ ipn{Xn)) (V?l) X ■ ■ ■ X . (4.1) 

Generally, this property breaks down when, in a statistical mechanical sense, the theory is in 
a mixed phase. In field theory language, this means there is more than one possible vacuum 
state. The clustering property is then restored by restricting the theory to the Hilbert space 
built on one of the vacua. In this sense, clustering is violated in a mild way, and to distinguish 
this from some other, potentially more serious, violations uncovered below, we will say that 
the theory satisfies a "generalized notion of clustering" . 

Let us consider the calculation of the Qn, the n-point function of the composite operator 
tr^A^: 

QniXi, . . . ,Xn) = {tT^\'^{Xi) ■ ■ ■tT^\'^{Xn)) . (4.2) 

For present purposes we restrict our attention to pure Af = 1 SU{N) gauge theory. Since 
tr^ A^ is the lowest component of a gauge-invariant chiral superfield (namely tr^ W"^ where 
W"' is the field-strength superfield), a well-known identity — reviewed in the Appendix — says 
that 

Qn{xi, . . . ,Xn) = const. , (4.3) 

independent of the Xj. Next let us consider this constant correlator in the instanton ap- 
proximation. This means that, at topological level k, X{x) is simply to be replaced by a 
general superposition of adjoint fermion zero modes in the general ADHM /c-instanton back- 
ground, weighted by Grassmann- valued parameters (i.e., fermionic collective coordinates). 
All bosonic and fermionic collective coordinates are then integrated over, in the appropriate 
supersymmetric way reviewed below. It can also be shown that Qn should still be a constant. 
(The field theory proof of the constancy of the correlation functions and its extension to the 
instanton approximation is discussed in the Appendix.) Now, in SU{N) gauge theory, at the 
topological level k, a multi-instanton has precisely 2kN adjoint fermion zero modes which 
need to be integrated over. Let us summarize the rules for Grassmann integration: if ^ is a 
Grassmann parameter, then 

jd^^ = l, jd^l = 0. (4.4) 
^For a discussion of clustering and other references, see Bogolubov et al. [Q. 
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Since tr^A^ is a Grassmann bilinear, it follows that Qn is only non- vanishing for n = kN. 
In particular, the one-point function Qi always vanishes. In summary, in the instanton 
approximation, at topological level k, we have the following selection rule: 

(tr^A2(x)) = Qi = for all A;; (4.5a) 



fc-inst 



fc-inst 



fc-inst 



(tr^A2(a;i) • ■ ■tr^A2(a;„)) = Gn , ^0 if and only if n = fcA^ . (4.5b) 



fc-inst 



Notice that these results already indicate a breakdown of clustering for the correlation func- 



tions (4^), although, as we shall explain below the breakdown is of the 'mild' variety and can 
be traced to the fact that in instanton approximation the theory is in a mixed phase, i.e. the 
instanton approximation samples the theory in a number of distinct vacua as opposed to a 
single vacuum. 



A general field-theoretic understanding of the selection rule ( [4.5a| )-( ^3E| ) was suggested 



in Refs. [0,^. The suggestion relies on the fact that, in A/" = 1 SU{N) gauge theory, the 
vacua of the theory come in an A^-tuplet [^. The vacua spontaneously break the discrete 
anomaly- free remnant of the classical f/(l)i? symmetry to the Z2 subgroup: Aq, — > — Aq,. 
The vacuum sector therefore consists of: 

{ \u) : < M < A^ - 1 } . (4.6) 
If we define the condensate J7 via 

J = (u = 0|tr^A2|M = 0) , (4.7) 
then the A^-tuple of vacua are related by phase factors, namely the A^**^ roots of unity: 

{u\tTj^X^\u) = J^e^™/^ . (4.8) 
Now let us see how the selection rule (^4.5a|) comes about. We define the density matrix 

Since the instanton calculation is 'Ij2n symmetric, it must average over all the vacua. This 
means 

N-l N-1 

(tr^ A^> = Tr(,tr^A^) = ^ E ^'n^' I") = ^ E = ^ ' (^.lO) 

Here the capitalized 'Tr' means a trace over the Hilbert space. In order to check the selection 
rule ( |4.5b| ), we need the additional assumption of a well-defined clustering limit. We have: 
(tr^A^(a;i) ■ ■ ■trjvA^(a;„)) = Tr(^tr^ A^(xi) ■ ■ ■ tr^ A^(a;„)) 

1 , , (4-11) 

= ^2Z<"I*^^^ (^1)^*^^^ (^2)P---Ptr^A2(x„)|M) , 



u=0 
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where P denotes the sum over a complete set of states. At this point the generahzed notion 
of the clustering assumption enters. We assume there exists a mass gap /i that is dynamically 
generated in the theory, and we consider the n insertion points are sufficiently far separated in 
Euclidean space compared to this scale: \xi—Xj \ ^ (Since Qn is a constant even in leading 
semi classical order, moving to this regime does not entail any additional approximations.) 
In this regime, the generalized cluster decomposition (in our present usage) is equivalent to 
the statement that P collapses to Pq where Pq is the projection operator onto vacuum states 
only: 

N-l 

P-.Po , Po = . (4.12) 

u=0 

Using the fact that the operator tr^ is diagonal in the u index,^ it follows that with the 
replacement (|4.12| ), the correlator (^4.11|) collapses to 



1 ^ \ 7^^ n - kN 

-V (Je^™/^)" = r (4.13) 
^ ^ ' 10 otherwise . 

Ji=0 1, 



Next we consider how this elementary analysis is modified if the iV-tuplet of vacua 
is supplemented by an extra vacuum state, the so-called Kovner-Shifman vacuum [0, which 
we denote \S). A single such vacuum is permissible under the discrete symmetry only if 

(^Itr^A^I^) = . (4.14) 

The analysis proceeds just as before, with the obvious modification that the density matrix 
Q should be replaced by g' , defined by 

N-l 

q' = {l-p)\S){S\ + , (4.15) 

where the probability p is a real number between and 1. Proceeding as before, we find that 
with the generalized clustering assumption 

N-l 

P^P'o, P[, = \S){S\ + , (4.16) 



n=0 



one derives 



(tr^A^(a:i)---tr^A^(a:„)) = <( ^^''^ (4.17) 

otherwise . 



^This follows from the fact that there should be no mixing between the sectors of Hilbert space built on 
each vacuum: they are super- selection sectors. 
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Obviously this modified expression also applies if there are several distinct KS vacua \Si) in 
which trA^ = 0, that is 



N-l 



d = + ' (4.18) 

where p = 1 — ^ O'l- 



N 

i=l u=0 



In the following we will calculate these (A;A^)-point correlators, first analytically for large 
N, then numerically for = 2 and k = 2, and will find a behavior quite different from either 
(i;Tl) or (j^). 



V Large- Calculation of Gluino Condensate Correla- 
tion Functions 

We now present an explicit evaluation of n = kN, in the limit N oo with k held fixed. 
Our answer turns out to be incompatible with both Eqs. ( [4.13| ) and (|4.17|) . The cleanest way 



to quantify this disagreement is to consider the (kNY^ root, (QkNY^''^ ■ In the large- limit, 
from Eq. (|4.13|) , i.e. clustering without the KS vacuum, one obtains: 

hm {GkNY^''' = J{N) . (5.1) 

N^oo 

We have written J' as J{N) to allow for an unknown dependence. Using the one instanton 
expression (|1.3| ) in the large- limit one expects 

lii^W-fAS (5.2) 

where e = 2.718 ■ • • . The key point is that the right-hand side of (|5.1| ) is independent of the 
topological number k (as well as of the space-time insertion points Xj). Note that Eq. ( p.l[ ) 
follows, not only from Eq. ( [4.13|) , but also from Eq. (|4.17| ), so long as the constant p either has 



a nonzero large- A^ limit, or else vanishes at large A^ more slowly than exponentially. Alterna- 
tively, with the "Shifman assumption" (|1.5| ) for p, which vanishes faster than exponentially 



at large A^, one obtains instead from Eq. ([4.17 



JtojW- = (|)"V(iV). (5.3) 



Combining this with the large- A^ limit of the 1-instanton expression (|1.3|) , one extracts instead 
the expression 

^ J(Ar) = A3 (5.4) 
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which now agrees (by construction) with the 1-instanton WCI calculation. 



Below we will calculate QkN, to leading order in but for all instanton number k, and 
will obtain a markedly different behavior. Explicitly we will find: 

^ hm {g.NY^''' = ^kA' + 0{N~') . (5.5) 



167r2 iv^oo'"''''' 

Notice that for A; = 1 we obviously recover the results (|5Tl|)-([57^) (or (p.3|)-(p^); however 
the linear k dependence is in sharp disagreement with the k dependence of either Eq. ( ^.1|) 
or Eq. ( |5.^j| ). This disagreement means that the generalized clustering assumption ( [4. 121 ) 



is invalid when combined with the instanton approximation. It also means that that the 
extension ( [4.16|) of this clustering assumption, in the presence of an extra KS vacuum state, 



is likewise invalid. 

The large- iV calculation proceeds as follows. | In supersymmetric theories, at topological 
level k, the bosonic and fermionic collective coordinates live, respectively, in complex-valued 
matrices a and A4, with elements: 

The indices run over 

M=l,...,Ar, ^,J = l,...,fc, d,/3 = l,2; (5.7) 

traces over these indices are denoted 'tr^v', 'tr^', and 'tr2', respectively. The elements of 
are Grassmann (i.e., anticommuting) quantities. The k x k submatrices a^^ = a'^<jj!ja and 
A^^ are subject to the Hermiticity conditions 

< = M'^ = M'^. (5.8) 

In the instanton approximation, the Feynman path integral is replaced by a finite-dimensional 
integration over the degrees of freedom in a and Ai. These /c-instanton collective coordinates 
are weighted according to the integration measure 1 11, 12,26, 2^0 

J ^phy^ \o\u{k) J ^ ^^g^ 

n \ n 5(itr,.r'^(tr2 r^aa)) J] 5 {iiuT'{Ma^ + d^M)Y 



X 

r=l,...,A;2 "c=l,2,3 a=l,2 



^Our conventions are taken from |T|,|T^ which also provide self-contained reviews of the ADHM formalism 
for the SU{N) gauge group. 

^°The reason we have 2^ rather than 2^'^' as in is that we restore Wess and Bagger integration 
conventions for the Ai' integration: J d?^ = 1 rather than 2 where = is the square of a Grassmann 
Weyl spinor. 
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where the two 5-functions enforce the bosonic and fermionic ADHM constraint conditions, 
respectively. The integrals over the k x k matrices a'^ and Ai' are defined as the integral over 
the components with respect to a Hermitian basis of k x k matrices T*" normalized so that 
tVkT^T'^ = 6^^. These matrices also provide explicit definitions of the 5-function factors in 
the way indicated. 



The form of the measure given in Eq. (|5.9| ) is known as the "fiat measure", since the 
bosonic and fermionic ADHM collective coordinates are integrated over as Cartesian vari- 



ables, subject to the nonlinear (5-function constraints. It was uniquely constructed in Ref. ^6 
to obey several important consistency requirements — including cluster decomposition — so 
that the failure of cluster uncovered below cannot be attributed to the collective coordinate 
measure. In practical applications, however, the fiat measure is not the most useful form 
available. When 

N > 2k , (5.10) 

it is convenient to switch to the so-called "gauge-invariant measure," involving a new set of 
variables in terms of which the arguments of the 5-functions are linear (and hence trivially 
implemented) [l^ . This is the form of the measure which we will utilize in the present section. 
The restriction ( p.lO[ ) is obviously well suited to the large- limit. As the name implies, the 



gauge-invariant measure can only be used to integrate gauge-invariant quantities, such as our 
present focus on correlators formed from tr^A^. Alternatively, for the special cases /c < 2, it 
is easy to solve the nonlinear constraints explicitly without such a change of variables P,p8 



In order to switch from the fiat measure to the gauge-invariant measure, one trades the 
collective coordinates w and w (which transform in the A^ of SU{N)) for the gauge-invariant 
bosonic bilinear quantity W, defined by |12| 



(l^"^)^^. = <«;„^.^, W' = tT2W, W' = ii2T'W, c=l,2,3. (5.11) 
The appropriate Jacobian for this change of variables reads: 

^2kN-^2kN^ _^ Cfc,^ (detsfciy)"^" W° J] d^"W\ (5.12) 



c=l,2,3 

where 

22kN-4k^+k ^2kN-2k^+k 

Note that the bosonic (5-function in ( [5.9| ) can be rewritten in a gauge-invariant way as the 
condition 



'^k,N = 2k — t; • (5.13) 



= W'^+K,a'J tr^ r^a"- = - , a'Jr^'^^ (5.14) 
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in terms of the gauge- invariant coordinates (here f/,J^^ is an 't Hooft tensor |29[). As advertised, 
these constraints are linear in the new variables W^] consequently the integrals simply 
remove the bosonic ADHM 5-functions in Eq. ( p.9[ ) (giving rise to the numerical factor of 
23fc2 fj,Qj^ |-]^g i'g ^]^g arguments of the 5-functions) . 

Next we perform a similar change of variables for the fermions, letting |T^ ] 

l^ui — U)ujce{C )ji ~l~ ^uij fJ'iu — iCa)ij^ju ~^ ^iu ' (5.15) 

where z/ hes in the orthogonal subspace to w: 

Wt^uj = , ^iuW^ja = ■ (5.16) 

One finds: 



j d^^fx d^^Ji W JJ 5 (trfc T"(A^a<i + d^M)) 



(5.17) 

where the (5-functions have been used to eliminate the C, variables from the problem. In 
summary, the gauge-invariant measure is: 

Vol [/(A;) 

and the constraint 5- functions have been eliminated for all k satisfying Eq. ( ^.10|) . For = 1, 



j d^'^'a' d'"w' d^'^'M' d^^\ d'''-''"ud'''''^'"u (det2fcl^)^-''= (5.1^ 



one recovers the expression ( p.l|) , a comparison which fixes the normalization constant Ci: 

C, = _2^+i/2a3^ . (5.19) 



For = 2, we recapture the Osborn measure discussed in Refs. [^^,|2^, which we utilize in 
Sec. VI below. 

Into this measure we now insert 

tr^A2(xi) X •■■ X tr^A2(xfciv) , (5.20) 

where the gluino A"(x) is replaced in the instanton approximation by a general superposition 
of adjoint fermion zero modes. In terms of the previously introduced collective coordinates 



a and Ai, a useful identity states [|l2 



tT^\\x) = -\ntT,MiV + l)Mf , (5.21) 
where the ADHM quantities V and / are defined as: 

P = l-A/A, f = {AA)-\ A = a + bx, (5.22) 
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and b is the (A^ + 2k) x (2A;) matrix whose lower 2k x 2k part is the identity and whose 
upper N X 2k part is zero (quaternionic muhiphcation is imphed in the product bx). As 
discussed earher, QkNi^i, . . . ,XkN) is actually a constant, independent of the Xj. The Xt can 
therefore be chosen for maximum simplicity of the algebra. However, the simplest conceivable 
choice, Xi = for all i, results in an ill-defined answer of the form "0 x oo" (the zero coming 
from unsaturated Grassmann integrations, and the infinity from divergences in the bosonic 
integrations due to coincident poles); we have already noted this fact in the 1-instanton sector 
in Sec. II above. The simplest choice of the Xi that avoids this problem turns out to be: 

Xi = ■ ■ ■ = XkN-k2 = , 

(5.23) 

XkN-k'2+l = ■ ■ ■ = XkN = X 

which we adopt for the remainder of this section.]^ 

In the large- iV limit the large preponderance of the insertions (|5.23|) are at Xi = 0, and the 
resulting factor of (tr^ A^(0))'^^~'^^ taken together with the Jacobian factor {det2kW)^~'^'^ 
from the measure ( |5.18| ), dominate the integral and can be treated in saddle-point approxi- 



mation. Below we will carry out this saddle-point evaluation in full detail, but we can already 
quite easily understand the source of the linear dependence on k in the final result ( |5.5| ) . The 
chain of argument goes as follows: 

1. Let us imagine carrying out all the Grassmann integrations in the problem. The 
remaining large- integrand will then have the form exp ( — A^F + (9(log A^)) where T might 
be termed the "effective large- A^ bosonic instanton action." The large- A^ saddle-points are 
then the stationary points of F with respect to the bosonic collective coordinates. By Lorentz 
symmetry, F can only depend on the four kx k matrices a'^ through even powers of a'^. (This 
is because the bulk of the insertions have been chosen to be at Xj = 0; otherwise one could 
form the Lorentz scalar x^a'^ and so have odd powers of a^.) It follows that the ansatz 

a^ = 0, n = l,2,3,4, (5.24a) 
^ = 0, c=l,2,3 (5.24b) 

is automatically a stationary point of F with respect to these collective coordinates. (Note 
that ( p. 24 bp follows automatically from ( [5.24a| ) by virtue of the ADHM constraints ( p.l4| ).) 



It will actually turn out that, once one assumes these saddle-point values, F is independent 
of the remaining collective coordinate matrix W^; furthermore we will verify that this saddle- 
point is actually a minimum of the Euclidean action. 
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As a nontrivial check on our algebra, we have also numerically integrated the large-iV correlator for 



insertions other than Eq. (5.23), and verified the constancy of the answer presented below. 
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2. Having anticipated the saddle-point ( [5.24a[ )-( p.24b| ) using these elementary symmetry 



considerations, let us back up to a stage in the analysis prior to the Grassmann integration, 
and proceed a little more carefully. Evaluating the insertions tr^A^(a;j) on this saddle-point, 
one easily verifies that the ( modes vanish when Xi = 0; consequently the ( integrations 
must be saturated entirely from the k"^ insertions at Xi = x. This leaves the v and v 
integrations to be saturated purely from the insertions at Xj = 0. Moreover, because M.' 
carries a Weyl spinor index a whereas v and v do not, the trjyA^(O) insertions depend on 
these Grassmann coordinates only through bilinears of the form u x u ot Ai' x Ai'] there are 
no cross terms. 

3. Performing all the Grassmann integrations then automatically generates a combina- 
toric factor 

{k')\{k')\{kN-2e)\ (^^^7^') . (5.25) 

Here the first three factors account for the indistinguishable bilinear insertions of the C, A^', 
and {v, z/} modes, respectively, while the final factor counts the ways of selecting the k"^ 
bilinears in Ai' from the kN — k"^ insertions at Xj = 0. Multiplying these combinatoric factors 
together, as well as the normalization constants Ck,N{Ci)'^ from Eq. (|5.18| ), and taking the 
{kNy^ root yields, in the large- limit: 

lim \ck,NiCif {k^. {kN - e)^'^''^ = 2^Ti^eN~^kk'^ + 0{N-^) . (5.26) 

Remarkably, apart from a factor of four, this back-of-the-envelope analysis precisely accounts 
for the previously announced final answer, Eq. ( |5.5| ). Note that most of the remaining 
contributions to the saddle-point analysis, which involve a specific convergent bosonic integral 
derived below, as well as the factor 2^^^^^ /Vol U{k) from Eq. ( |5.18D , reduce to unity when the 
{kNy^ root is taken in the large-iV limit; the missing factor of four will simply come from 
the leading saddle-point evaluation of the bosonic integrand. 

Here are the details of the large- calculation of QkN- Since the problem has an obvious 
U{k) symmetry [T^, we will find it convenient to work in a basis where (which transforms 
in the adjoint of the U{k)) is diagonal: 



/2p? \ 



V 2plJ 



(5.27) 



As the notation implies, in the dilute instanton gas limit pi can be identified with the scale 
size of the i^^ instanton in the fc-instanton sector (see Sec. II. 4 of |jl2|). The appropriate 
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change of variables reads: 



Vol U{k) 



dpi- -dpi n (P'-P?)'- (5-2^ 



For k = 1 one has, of course, J dW'^ — > 2 dp^. 



Now let us consider the Grassmann integrations, beginning with the ( modes. We assume 
the saddle-point conditions ( ^■24a| )-( |5.24b ), in which case 



A 



w 



X ■ 1 



[felx[fel 



and from Eq. ( p.21D 



/ 1 



f 



\ 



Pk 



(5.29) 



tr^A^(a;) 



where 



ID ^ 



(5.30) 



(5.31) 



and the omitted terms in Eq. (|5.3CI|) represent dependence on the other Grassmann modes 
{Ai', z/, u}. It is obvious from Eq. ( [5.31| ) that -Fij(O) = 0, so that the ( modes must be entirely 
saturated from the k"^ insertions at claimed above. Performing the ( integrations 

then yields 



-ifik'y. l[F,,{x). 



(5.32) 



Next we consider the insertions at Xj = 0. Focusing on the A4' modes first, one finds from 
Eq. dm): 



tr^A^(O) = 2j2iM'%AM'^hip;' + Pj' + pfpf) + 
omitting the z/ x z/ terms. Hence the Ai' integrations yield 

k 



(5.33) 



kN-k'^ 
k^ 



(5.34) 
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where the combinatoric factors in ( |5.34| ) (as well as in ( |5.32|) ) have been explained previ- 
ously.0 



Finally we turn to the {z/, u} integrations. Since (unlike the ( and Ai' modes) the num- 
ber of u and u modes grows with as kN — it does not suffice merely to plug in the 
saddle-point values ( p.24a )-( 5.24b ) and ( ^.291 ). One must also calculate the Gaussian deter- 
minant about the saddle-point, which provides an 0{N^) multiplicative contribution to the 
answer. Accordingly we expand about ( 5.24a|) - (|5.24b| ) to quadratic order in the a'^. The uxu 
contribution to tr^A^(O) has the form 



2 u ^ui n fij 



x=0 



x=0 



(5.35) 



as follows from Eqs. ( p.21D -( [5.22| ), and Eq. (2.63) of Performing the integrations 
therefore gives 

{kN - 2ky. exp ({N - 2k)tT. log (2/ ■ tr, bVb ■ f) I 

V \x=0 

k 4 

(fcAr-2fc2)!exp((Ar-2A;)(logdeta6(iy°)-2 - | (P."' + Pj') + OKY) 

i,j=l n=l 

(5.36) 



The negative sign in front of the quadratic term in a'^ confirms that our saddle-point (^.24a|) - 
(|5.24b| ) is in fact a minimum of the action. Combining this expression with the measure 
factor in Eq. ( 5.18| ), namely 

(detsfcVr)"^-'^ = exp((Ar-2A;)logdet2fcVr) 

= e^p({N-2k){logdetk{lWy + 0«)^) 



(5.37) 



and performing the Gaussian integrations over a^, yields: 

k 
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where the omitted terms are suppressed by powers of A^. 



+ 



(5.3^ 



Finally one combines Eqs. (^^), ( F^ ), (1041) and ( ^^Hj) to obtain the leading- 

order result for the correlator: 



lim GkN 

N^oo 



25kN+k^-k+l/2 ^2kN~k+l/2 ^kN j^^2-ji f^kN~k'^+l/2 ^ 
32fc2 ]\fkN+k2-l/2 



3kN 



(5.39) 



^^One can easily check that these large- iV formulae are consistent with the explicit 1-instanto n calc ul ation 
presented in Sec. II which is exact in N. In particular, if one takes Xi = x while Xj = 0, then Eqs. (2.5b)-(2.5c) 



are suppressed vis-a-vis Eq. ( 2.5a ) by factors of a^, and in turn, ^ N as follows from Eq. ( 5.36 ) below 
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where is the convergent integral 



= dpi-., dpi n iPl-p]f■\{F^A^){^-{P^/P^+P^IP^r')■ (5.40) 

l<i<j<k i,j=l 



Note that is independent of x as a simple rescaling argument confirms. For the simple case 
k = 1, the {pf — pjY terms in this integral are absent; one finds Xi = | and the expression 
( |5.39|) agrees — as it must — with the large- limit of the 1-instanton SCI result (|2.11| ). 



VI The 4-Point Function of the Gluino Condensate in 
SU(2) Gauge Theory 

We have seen that cluster decomposition fails (both with and without a KS vacuum) in the 
SCI calculation of the gluino condensate, for gauge group SU{N) in the large-A^ limit. In 
this section we focus instead on the gauge group SU (2). In this case, at the 1-instanton level, 
the 2-point function (|1.2| ) works out to: 



trA2(xi) trA2(x2)\ 
167r2 167r2 / 



fA^ (6.1) 



Here we will calculate the 4-point function, which receives a nonzero contribution at the 
2-instanton level: 



cA^' . (6.2) 



trX'^ixi) trA^(x2) trX'^ix^) trA^(x4) 

167r2 16^^2 ^6^2 IQ^2 

In the absence of a KS vacuum, generalized cluster decomposition together with Eq. ( |6.1|) 
predicts c = (4/5)^ = .64. Alternatively, in the presence of a KS vacuum, weighted according 
to Eq. (|1.5|) in order to reconcile the SCI and WCI 1-instanton calculations, one expects 
c = 4/5 = .8. Instead, we have calculated c numerically, and find: 

c ~ .500 ± .026 . (6.3) 

Here are the details of the calculation. 

As mentioned above, for k = 2, one can eliminate the (5-function constraints in Eq. ( p.9| ) 
without changing variables. Another simplification for the particular gauge group SU (2) is 
that one can adopt a concise quaternionic representation for the ADHM bosonic collective 
coordinates, taking advantage of the fact that SU{2) = Sp{l). Specifically, the 16 gauge and 
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8 gaugino collective coordinates live, respectively, in the following matrices:^ 

M'u, M[,A , (6.4) 

where and the matrices a„ as well as A^-y are real- valued (unlike the 

complex-valued collective coordinates of the same name introduced in Eq. ( |5.6|) which are 
needed for general SU{N)). The resulting 2-instanton "Osborn measure" on these collective 
coordinates is detailed in Refs. [^|2^,^, and reads: 



j dfil^^^ = 2^^A^ j d^wi d*W2 d*a[^ d*a'^2 ^^Vi d^l^2 d^M[^ d^M'^^ ^ |2^^ ' ^^'^^ 



Here the (5-function constraints from the fiat measure have been used to eliminate ^^^1 
A^'i9 in terms of the other collective coordinates, via: 



a' 



12 ^ ^^4(^2^1 -wiwa) , (6.6) 



and 



•^'12 = T^TTyi 0,3 ( Sa'is-^S + «^2/Ul - Wi/i2 ) , (6.7) 



where we have defined 



M', = HM'u-M'22) . (6.8) 



Into this measure one inserts the 4-point function of the classical condensate, expressed as 
a function of the 2-instanton collective coordinates ( p.4|) . The 8-dimensional Grassmann 
integrations over {/ii, /i2, Ai'i^, 7^22) then accomplished in two steps. The first step 
is to expand the integrand in terms of Grassmann variables using a modified version of 
the program "Dill", written for Mathematica.0 The second step involves the explicit 
Grassmann integration, accomplished using an "awk-script" implemented on a UNIX system 
and made to perform the symbolic algebra of Grassmann integration. 

The resulting 16-dimensional bosonic integration over {wi, W2, a'n, 022)5 remaining 
quaternionic variables, is carried out using a standard Monte Carlo integration procedure. 
The integrable singularities are handled using the standard procedure: firstly, dropping a tiny 
region around the integrable singularities and then making sure that the contribution from 



^•^See Ref. for details of notation and conventions pertinent to Sec. VI. 

14 "j3j^" jg Mathematica package originally written by Vladan Lucie in 1994 in order to simplify 
SUSY algebraic expressions. This program can be modified so that it can handle the large number of 
Grassmann variables that we need. 
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this dropped region is negligibly smaller than the precision required. After 450 million points 
have been sampled, we have obtained the numerical value (|6.3|) given above. As a check 
on our numerics, we have also verified the constancy of the answer by comparing different 
choices for the four space-time insertion points. 



VII Discussion 



The mismatch between the strong coupling and weak coupling calculations is a fascinating 
puzzle. Previously, only the mismatch at the one instanton level was known; now we see a 
mismatch established at large for all instanton numbers, and for A^ = 2 at the 2-instanton 
level. Certainly we do not mean to imply that, because of this mismatch, SCI calculations 
are all necessarily suspect; indeed an A/" = 4 supersymmetric version of an SCI calculation 



performed by some of us ||12| has recently provided a dramatic quantitative and qualitative 
verification of Maldacena's conjecture. However in this case the coupling does not run and 
the calculation can be performed at weak coupling, actually small g'^N , where the instanton 
approximation is fully justified. The continuation to strong coupling, large g'^N , is then 
accomplished by means of a non-renormalization theorem. Rather, our objections to the SCI 
computation are more narrow and technical in scope: specifically, our calculations imply a 
fundamental breakdown of clustering in the instanton approximation to the gluino condensate 
at strong coupling. 

One may wonder what the origin for this breakdown is? The usual justification for the 
strong coupling calculation is that one can take \xi — Xj \ much smaller than the scale of strong 
coupling effects and so the theory would be weakly coupled, due to asymptotic freedom, 
and the instanton calculation would be justified. Then, since the correlation functions ( [4.2|) 



are independent of the positions, the result would be valid at all distances. This point-of- 
view has simultaneously been used and criticized by various authors The asymptotic 

freedom argument means that the first-order, second-order, etc., perturbative corrections 
to the SCI calculations are small — indeed, for the gluino condensate correlators discussed 
herein, these perturbative corrections are entirely absent due to a nonrenormalization theorem 
[0. However, the asymptotic freedom argument does not guarantee that the zeroth order 
instanton calculation is itself complete; there may be other non-perturbative configurations 
contributing to the correlators. Indeed, the breakdown of cluster suggests that such additional 
nonperturbative configurations (with size of order A~^) must be present, and that they must 
account for the mismatch between the SCI and WCI calculations. 
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In contrast, it seems that the WCI calculation uses a method that has amassed a con- 
siderable pedigree. These kinds of calculations appear to be consistent in all applications 



and agree with other non-instanton methods pO[; for example, the two-instanton check of 



the Seiberg-Witten approach to A/" = 2 theories P,|T3] and the latter calculation in Sec. III. 
Moreover, in the WCI set-up, large-scale nonperturbative configurations as just discussed, 
such as instanton-antiinstanton pairs of size A~^, would be exponentially suppressed in the 
path integral so long as A ^ v. We should further note that, as the separation between 
insertions tends to zero, the WCI calculation does not smoothly go over to the SCI calcula- 
tion as one would naively expect; there are additional important contributions which will be 
discussed in a separate publication (work in progress). 

It is unfortunate that, based on our results, the highly original and intriguing (both 
theoretically and phenomenologically) proposal of Shifman, namely the existence of a chirally 
symmetric vacuum state, loses much of its raison d'etre. Then again, we have not actually 
ruled out the existence of such a state. After the completion of this work, it has been 
suggested that the mixing parameter p of the KS vacuum, defined in Eq. ( [4.15|) above, may 
actually be instanton number dependent Prima facie, this appears to be incompatible 



with invariance under large gauge transformations (|/c) |A; + 1)); however, if such a counter- 
intuitive flexibility is permissible in the definition of the instanton vacuum, clustering in the 
presence of the KS vacuum can be saved. 

Conceptual difficulties with the instanton approximation and cluster decomposition were 
pointed out in the context of pure (non-supersymmetric) QCD some time ago. Since this may 
have some bearing on the present discussion, we review some comments of Liischer regarding 



this issue [^. The pure instanton (i.e. no anti-instantons) approximation to QCD obviously 



violates parity since 

(trF„„ = p 7^ , (7.1) 

where p here is the instanton density. Parity is then recovered by summing over instantons 
(I) and anti-instantons (I); however, in this approximation the cluster property would not 
hold. To see this note that 

(trF„„*F„„(x) trFp/Fp,(0)) p V , (7.2) 

\x\—*oo 

whereas 

In order to resolve this clustering conundrum, it is apparent that additional configurations, 
which may, in the dilute gas limit, be thought of as mixtures of instantons and anti-instantons, 
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would need to be incorporated in the approximation. In this case the two-point function ( [7.2[ ) 
would indeed be zero, the result of summing the II, II, II and II contributions which on 
average all contribute equally. Away from the dilute instanton gas limit, the identification and 
physical interpretation of these additional cluster-restoring nonperturbative configurations is 
necessarily more subtle. 

In summary, the results of this paper imply something analogous in the = 1 theory: 
additional configurations must contribute to the correlators at strong coupling and resolve 
the breakdown of clustering (as well as repairing the mismatch between the SCI and WCI 



calculations). In fact, it was suspected some time ago (see Ref. p8|,^ and references therein) 
that in strongly coupled theories, it may be more appropriate to think of instantons as com- 
posite configurations of some more basic objects: so-called "instanton partons". The domi- 
nant contributions to the path integral at strong coupling would then arise from the partons 
themselves. In Ref. , we make this piece of folklore more precise by identifying instanton 



partons with the monopole configurations of the supersymmetric Yang-Mills compactified on 
the cylinder x 5*^, with the circle having circumference (3. Each monopole has precisely two 
gluino zero modes, rather than four for the instanton. The instanton itself is then identified 
with a specific two-monopole configuration. We calculate the monopole contribution to the 
gluino condensate and then, at the end of the day, take the decompactification limit (3 ^ oo. 
The value of the gluino condensate obtained in this way, is precisely the WCI result (|1.4|). 
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Appendix A: Constancy of the Correlators 



Since the constancy of the correlation function ([4.2|) plays such an important role in our 
analysis, in this appendix, we review the arguments leading to this result. More importantly, 
we explain how the field theory proof remains valid in the instanton approximation. 
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First the field theory proof [|T],0]. The argument is completely general and applies to the 
correlation functions of any lowest component A of a gauge invariant chiral superfield 

^ = A{x) + V2eilj{x) + ■ ■ ■ . (A.l) 

In the present discussion, the operator tr is the lowest component of the chiral superfield 
triy^PFo;, where Wa is supersymmetric field strength. Consider the correlation function 

(Aiixi) ■ ■ ■ Ap{xp)) . (A.2) 

We will show that this is independent of the Xj's. To this end, one has 



Hence 
d 

Q-^ (M^i) ■ ■ ■ M^p)) 

=|<"(0|Ai(a;i) ■ ■ ■ Ai_i(xi_i){Q<i, ^pia{x^)}Ai+l{x,+l) ■ ■ ■ Ap{xp)\0) 

= - i^n" ■ ■ ■ [^°' ^■(^i)] ■ ■ ■ A-lix^-l)^^aiXi)Ai+l{x^+l) ■ ■ ■ Ap{Xp)\0) 

P 

(0|^i(a;i)---A,_i(x,_i)7/^i«(xi)Ai+i(xi+i)---[QA,Aj-(xj)]---Ap(xp)|0) , 

j=i+i 

(A.4) 

where the last line follows by commuting the Qa through the other insertions, to the left and 
right, respectively, until it hits the vacuum which it annihilates. But [Qo, = and 

therefore the right-hand side of ([A.4|) vanishes and consequently the correlation function is. 



indeed, independent of the insertion points. QED. 

If the multiple correlator of tr is constant in the full field theory, it then becomes an 
issue as to whether this constancy is retained in the instanton approximation. That it is, 
rests upon two facts. Firstly, the supersymmetry transformations of the fields can be traded 
for supersymmetry transformations of the collective coordinates 10, 11, 26, ^|. In other 



words, the supersymmetry algebra is represented on the collective coordinates. Specifically, 
under an infinitesimal supersymmetry transformation + S,Q'- 

Saa = laM, = -MC (A.5a) 

5M = -4ib"U = AiCK ■ (A. 5b) 



In particular, ( |A.3[ ) will hold, with the fields replaced by their expression in the instanton 



background and with the right-hand side involving the appropriate transformation of the 
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collective coordinates. Ipso facto, the argument leading to ( |A.4[ ) will hold with the transfor- 
mations acting on the collective coordinates; moreover [Qai ^] = 0, understood as a trans- 
formation of the collective coordinates. The remaining piece of the proof is the analogue 
of the fact that Qa annihilates the vacuum state. In the instanton approximation, where 
the functional integral is approximated by the integral over the collective coordinates, the 
analogue of the statement that the vacuum is a supersymmetry invariant, is the statement 
that the measure on the space of collective coordinates is invariant under the supersymmetry 
transformations (|A.5a|)-([A.5b|). This invariance was proved in pGlp^. 
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